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Abstract 

Perceptual learning has been used to probe the mechanisms of cortical plasticity in the adult brain. Feedback projections 
are ubiquitous in the cortex, but little Is known about their role in cortical plasticity. Here we explore the hypothesis that 
learning visual orientation discrimination involves learning-dependent plasticity of top-down feedback inputs from higher 
cortical areas, serving a different function from plasticity due to changes in recurrent connections within a cortical area. In a 
Hodgkin-Huxley-based spiking neural network model of visual cortex, we show that modulation of feedback Inputs to VI 
from higher cortical areas results in shunting Inhibition in VI neurons, which changes the response properties of VI 
neurons. The orientation selectivity of VI neurons is enhanced without changing orientation preference, preserving the 
topographic organizations in VI. These results provide new Insights to the mechanisms of plasticity in the adult brain, 
reconciling apparently Inconsistent experiments and providing a new hypothesis for a functional role of the feedback 
connections. 
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Introduction 

The adult brain remains plastic long after the developmental 
period [1]. This ability to remain plastic is fundamental for the adult 
brain to be able to learn and adapt in the ever-changing sensory 
environment. However, we do not fuUy understand the limits and 
mechanisms of the adult brain plasticity [2] . Perceptual learning - 
improvement in perception due to experience with stimuli - has 
been used to explore plasticity in sensory cortices [3] . In particular, 
stimulus specificity of the improvements of orientation discrimina- 
tion [4—6] suggests that the primary visual cortex, area VI, where 
neurons have small receptive fields, is the site for cortical changes 
that underlie the improvement of orientation discrimination. 
However, it is still unclear what cellular and synaptic changes in 
VI are responsible for such perceptual learning. 

Reports of learning-dependent changes of VI neurons are 
inconsistent: Learning to discriminate orientations of visual stimuli 
in one study resulted in sharpening orientation tuning curves in a 
subgroup of VI neurons [7] (Fig. lA-D); others found instead that 
responses to the trained orientation reduced in VI neurons [8] 
(Fig. 1E,F). In both studies, VI neurons responded preferably to 
the same orientations as they did before learning, thus preserving 
the topographic organization in VI. 

Previous models of visual perceptual learning based on plasticity of 
the recurrent and feedforward connections [9-10] could neither 



explain the observed stability of orientation preferences of VI neurons 
nor reconcile the fmdings of the aforementioned experiments. 

Here we suggest that learning-dependent changes in VI cortex 
involve top-down projections into VI from higher cortical areas. 
In a spiking neural network model of visual cortex (Fig. 2), we 
showed that repeated stimulus presentations resulted in strength- 
ening the feedforward connections from VI neurons to neurons in 
a higher visual area (such as area V2). Greater activity in V2 in 
turn led to strengthening the feedback connections from V2 to VI 
neurons, which helped to maintain the stability of the V1-V2 
network through shunting inhibition [11-12]. As a consequence of 
the shunting inhibition, VI neurons had enhanced orientation 
selectivity and responded to a narrower range of orientations. The 
diffiise nature of the feedback inputs to VI neurons allowed 
improvement of orientation selectivity in VI neurons without 
changing their orientation preferences, thus preserving the map of 
orientation representations in V 1 . 

The feedback model reproduced both the previously observed 
sharpening and reduction of neural responses, but not in the same 
VI neurons. We assume that previous studies might have recorded 
from different populations of VI neurons. We suggest that the 
stability of cortical representations in VI may be critical for 
normal visual processing, and learning-dependent plasticity of 
visual cortex could be based on mechanisms that preserve the 
topography of cortical maps. 
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Author Summary 

The specificity of visual perceptual learning suggests the 
primary visual cortex, area VI , is the site of neural plasticity 
that underlies learning. However, VI does not exhibit 
topographic changes observed in other sensory cortices 
with learning, and changes observed in VI following 
perceptual learning are contradictory. In a spiking neuro- 
nal model of cortical network, we show that feedback 
inputs to VI can change the response properties of VI 
neurons without changing the topographic organization in 
VI. The feedback model reconciles previous experiments 
and provides a new hypothesis for a functional role of the 
feedback connections. 

Results 

Shunting inhibition of VI neurons 

To explore mechanisms of leaming-induceci changes in V I we 
designed a simple network model consisting of two cortical areas 
(VI and V2) and including plastic feedforward and feedback 
connections (Fig. 2; see also Models). In the model, a stimulus was 
presented to VI neurons on many repetitions, which resulted in 
strengthening the feedforward connections from V 1 to V2 neurons 
(Fig. 3A). Higher levels of activity in V2 neurons in turn led to 
strengthening the feedback connections from V2 back to VI 
neurons. The excitatory feedbacks targeted both the excitatory 
and inhibitory neurons in VI (Fig. 3B) and resulted in balanced 
recurrent excitatory and inhibitory currents in V 1 , which canceled 
each other at the resting potential of VI neurons leading to 
shunting inhibition of VI neurons (Fig. 3C). Strengthening the 
feedback inputs with learning did not change the balance but 
increased the amplitudes of both the recurrent excitatory and 
inhibitory currents (Fig. 3D). Thus, although the feedback inputs 
were excitatory, they suppressed activity in VI neurons through 
shunting inhibition mechanism, thus maintaining stability of the 
network. 

Due to the diffuse nature of the feedback inputs in the model, 
the shunting inhibition affected responses of all VI neurons. By 
applying a range of input stimuli, we tested f-I curve in a VI 
neuron that did not respond to the trained stimulus but received 
the modified feedback input. After training, the f-I curve of the V 1 
neuron shifted to the right, resulting in a subtractive effect on 
neural responses (Fig. 3E). This was because the strong feedback 
inputs resulted in stronger recurrent excitatory and inhibitory 
input currents, which decreased input resistance at the resting 
potential of the neuron. As we the depolarized membrane 
potential by applying a direct current to the VI neuron, the 
driving force of the inhibitory current increased and the driving 
force of the excitatory current to the neuron decreased causing 
shunting inhibition in the neuron. Thus, shunting inhibition can 
reduce neural responses by counteracting the driving stimulus, as 
was previously predicted [12]. 

Different enhancement of orientation selectivity in VI 
neurons after learning 

After training the model with a stimulus, we observed that both 
amplitude and width of tuning curve of a VI neuron that preferred 
the trained orientation reduced due to the shunting inhibition 
(Fig. 4A,B). However, slope of the neuron's tuning curve did not 
change significandy (Fig. 4D). We refer to this change in the 
tuning curve as "reduction" of the response. In another VI 
neuron, for which the trained orientation was not the preferred 
one, reduction of the width of the tuning curve was greater than 



reduction of the response amplitude (Fig. 4B,C), which steepened 
the slope of the tuning curve (Fig. 4D). We refer to this change as 
"sharpening" the tuning curve. 

Thus, the reduction of responses was observed in VI neurons 
that preferred the trained orientation, and sharpening of tuning 
curves was observed in VI neurons that responded to the trained 
orientation but preferred orientations different from the trained 
one. 

Stimulus specificity of tuning curve changes 

To understand causes of different changes of the tuning curves 
of VI neurons (Fig. 4), we explored stimulus dependence of neural 
changes in V2. Learning resulted in strengthening connections 
from excitatory neurons representing the trained stimulus in V 1 to 
excitatory neurons representing the trained stimulus in V2 
(Fig. 5A). Higher activity in those V2 neurons resulted in strong 
feedback inputs to all VI neurons (Fig. 5D), which led to overall 
reduction of neural responses in VI due to shunting inhibition 
(Fig. 5E, red lines). A novel stimulus activated a population of V2 
neurons that included some V2 neurons with previously modified 
connections and other V2 neurons with unmodified connections 
(Fig. 5B). Hence, the feedback inputs from the experienced V2 
neurons were combined with feedback inputs from naive V2 
neurons, which resulted in a lesser increase of the overall feedback 
inputs to VI neurons (Fig. 5D) and led to a sharper tuning curve 
(Fig. 5E, green lines). Finally, a stimulus that activated V2 neurons 
with unmodified synapses (Fig. 3C) did not change the feedback 
inputs to VI (Fig. 5D), therefore, no changes of tuning curves were 
observed in VI neurons (Fig. 5E, purple lines). Thus, after 
learning, strength of the feedback input depended on orientation 
of a stimulus: The feedback inputs were the strongest for the 
trained orientation and weakened as stimulus orientation moved 
away from the trained orientation. 

Effect of feedbacks on response functions in VI 

To understand how varying strength of the feedback inputs led 
to either firing rate reduction or sharpening of the tuning curves of 
VI neurons (Fig. 4A and 4C), we analyzed how recurrent 
interactions within VI network change response properties of 
VI neurons by comparing the responses of two VI neurons 
positioned at different locations with respect to the stimuli. 

In agreement with experimental data from VI [13-15], 
recurrent interactions among VI neurons increased responses of 
strongly responding VI neurons and suppressed responses of 
weakly responding VI neurons (Fig. 6A). This occurred because 
recurrent interactions effectively increased nonhnearity of input- 
output relations in VI neurons (Fig. 6B). The new input-output 
function had an intermediate area where the slope changed most 
strongly. W e refer to this region as an area of strong nonhnearity. 

The trained stimulus, which was the preferred stimulus for 
neuron A (Fig. 6C), evoked strong feedback inputs and resulted in 
proportional reduction of neural responses in both strongly (e.g. 
neuron A) and weakly (e.g. neuron B) responding VI neurons 
(Fig. 6C). This proportional reduction was because the strong 
feedback equally affected neurons A and B by shifting responses in 
both neurons beyond the area of strong nonhnearity (Fig. 6D). 

A novel stimulus, which was the preferred for neuron B, evoked 
weaker feedback inputs and led to the sharpening of VI neural 
responses (Fig. 6E). This sharpening was because not-so-strong 
feedback affected neurons A and B differentiy. The feedback 
shifted the response of neuron A beyond the area of strong 
nonhnearity, which led to a strong reduction of its response. In 
contrast, the response of neuron B remained outside the area of 
strong nonhnearity and displayed only a minor decrease (Fig. 6F). 
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Figure 1. Effect of training on orientation discrimination in V1 neurons. A. Orientation discrimination task. Subjects reported whether the 
test orientation was tilted clockwise or anticlockwise with respect to the trained orientation (from [7]). B. Performance in orientation discrimination 
task (from [7]). C. Orientation tuning curves of five sample VI neurons (from [7]). D. Slope measured at the trained orientation for trained neurons 
(solid red line) and naive neurons (dashed blue line) (from [7]). E. Performance in orientation discrimination task in another experiment (from [8]). F. 
Neural responses in VI in the trained location (grey line) and an untrained location (black) (from [8]). 
doi:10.1371/journal.pcbi.1003770.g001 



This difFerence in the firing rate change of different V I neurons 
resulted in sharpening of the population response (Fig. 6E). 

Thus, the strong nonUnearity of the input-output function of VI 
neurons and the varying strength of the feedback inputs produced 
different changes in the neural population, reducing the responses 
of some neurons and sharpening the tuning curves of others. These 
changes of the population responses observed in our model were 
similar to the reported changes of the tuning curves in V 1 neurons 
following perceptual learning [7-8] . 



Discussion 

In the primary auditory, somatosensory, and motor cortices, 
learning leads to dramatic increases in the number of neurons that 
represent the trained stimuli [16-20]. In the primary visual cortex 
(VI), topographic reorganization was observed following severe 
sensory changes [21-23], which raised the question of whether 
such changes in VI can also occur during normal perceptual 
learning [24]. Topographic reorganization in VI have not been 
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Figure 2. Model of visual cortex. Layer one represented layer 2/3 of 
VI and contained excitatory and inhibitory neurons. Layer two 
represented V2 area and contained excitatory neurons. VI neurons 
were connected through recurrent (excitatory and inhibitory) connec- 
tions. VI neurons received bell-shaped input, which mimicked 
orientation-tuned input from layer 4 neurons of VI. VI excitatory 
neurons sent convergent feedforward projections to excitatory V2 
neurons. V2 neurons projected back to all VI neurons. Stimuli of 
different orientations were modeled by shifting the center of the input 
along one-dimensional network. 
doi:10.1371/journal.pcbi.1003770.g002 

reported in any study following perceptual learning. Here we 
have shown cellular and synaptic mechanisms of cortical 
plasticity that can preserve the topographic organization in area 
VI while allowing changes to the properties of VI neurons. 
These results suggest that learning to discriminate simple visual 
stimuli may evoke changes that are different from the 
topographic reorganization observed in other sensory cortices. 
The feedback mechanism for learning proposed here may not, 
however, be limited to the visual system. Other sensory modalities 
may share similar mechanisms, depending on the particular 
learning task. 

Scope of the model 

Our model pertains to experiments that involve learning to 
discriminate low-level features of visual stimuli presented at a 
particular location within receptive field of VI neurons. We 
explored here if learning-induced changes in V 1 were not due to 
synaptic plasticity of short-range recurrent connections among 
nearby V 1 neurons within a hypercolumn but instead were due to 
changes in feedback inputs to VI neurons from neurons in higher 
visual areas. There may be other types of plasticity in the long- 
range horizontal connections between distant hypercolumns that 
can explain other types of learning experiments involving complex 
tasks and large stimuli [25-27]. 

Robustness of the model 

Reduction of neural responses and sharpening of tuning curves 
in our model resulted from strong nonlinearity of input-output 



relations in VI neurons and varying strengths of stimulus evoked 
feedback inputs to VI neurons. The intrinsic nonlinearity of input- 
output relations in VI neurons was further amplified by recurrent 
interactions. The variation of strengths of the feedback inputs 
resulted from the assumption that a population of V2 neurons 
representing a stimulus had a finite size. Therefore, the trained 
and novel stimuli evoked responses in an overlapping but different 
neuronal populations, which resulted in diflFerent feedback inputs 
to VI neurons upon presentation of the trained vs. novel stimuli. 
Thus, both reduction of the neural responses and sharpening were 
observed in VI network but in different populations of VI neurons 
with respect to the location of stimulation. The learning process in 
the model was incremental: More training resulted in stronger 
feedforward and feedback connections, and therefore led to 
stronger reduction and sharpening. Robustness of the results was 
tested by varying the synaptic weights in the model, which 
produced consistent behavior over a range of parameter values. 
Other aspects of the model were based on neural data and did not 
affect the conclusions. 

Shunting inhibition and gain modulation 

In a previous study, driving a spiking neuron with balanced 
excitation and inhibition currents changed the slope of the gain 
function (input-output function) [11]. The changes of the gain 
functions were caused by shunting inhibition due to balanced 
excitation and inhibition and noise in the input currents. However, 
in perceptual learning experiments analysis of noise in VI neural 
responses before and after learning did not show any change in the 
noise level [7]. In the current model, input noise did not change 
with training because the firing rates in VI did not increase with 
learning, and firing rates in V2 were maintained through the 
feedback mechanism. Thus, in the model, balanced excitatory and 
inhibitory feedback currents were solely responsible for the 
subtractive inhibition in VI neurons, in agreement with the 
experiments. 

Learning and adaptation 

In our model, the feedback mechanism of learning resulted in 
symmetric changes of the tuning curves in V 1 . However, tuning 
curves of V 1 neurons can in fact adapt and become asymmetric in 
response to repeated stimulus presentation [28]. This adaptation 
may be due to short-term synaptic depression and facilitation in 
the recurrent connections, which can result in asymmetric changes 
of the tuning curves. In contrast, the changes in VI discussed in 
our study are due to long-term changes in the feedforward and 
feedback connections. 

Both adaptation and learning could take place during stimula- 
tion and two processes could interact with each other [29-31]. 
The critical distinction between adaptation and learning is that 
learning, unlike adaptation, is task dependent. Both adaptation 
and learning can take various forms due to diversity of plasticity 
mechanisms (Hebbian, homeostatic, reinforcement, etc) that could 
be evoked by the stimulus [32]. However, with learning to perform 
a task, changes that enhance performance are reinforced, and 
changes that do not improve performance are suppressed, a 
distinction that may depend on a reinforcement signal. We have 
not included a reinforcement signal in our model, and instead 
assumed that plasticity of feedforward and feedback connections 
somehow depended on the task. Although recurrent connections 
did not contribute to learning in our model, they could be the basis 
for adaptation. Thus, both learning and adaptation can occur 
simultaneously and their relative contribution can be context 
specific. 
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Figure 3. Shunting inhibition of VI. A. With learning, feedforward connections from a population of VI neurons to a population V2 neurons 
strengthened. B. The excitatory V2 neurons projected to both the excitatory and inhibitory neurons in VI. C. Top. Membrane potential of a VI 
neuron before training. Bottom. Excitatory and inhibitory input currents to the VI neuron before training. D. Top. Membrane potential of the VI 
neuron after training. Bottom. Excitatory and inhibitory currents to the VI neuron after training. E. Response function (f-l curve) of a VI neuron. The 
neuron did not respond to the trained stimulus, and the response function was tested by injecting current to the neuron. 
doi:1 0.1 371 /journal.pcbi.1 003770.g003 
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Figure 4. Enhancement of stimulus selectivity. A. Tuning curve of a VI neuron that preferred the trained orientation; before (black) and after 
(green) training. B. Amplitude-to-width ratio (R) of the tuning curve before (black) and after training (green) in a neuron that preferred the trained 
orientation, and in a neuron that preferred different orientation (blue). C. Tuning curve of a VI neuron that preferred orientation different than the 
trained one; before (black) and after (blue) training. D. Slope change of the tuning curve of the VI neuron the preferred the trained orientation 
(green) and another VI neuron that preferred another orientation (blue line with filled circles). 
doi:1 0.1 371 /journal.pcbi.1 003770.g004 



Specificity of perceptual learning 

In many experiments on visual perceptual learning, the 
improvements of perception were stimulus specific [4-6], 
suggesting that early visual areas were possible sites for the 
changes underlying visual perceptual learning. However, 
under some conditions learning can be generalized to other 
locations in the visual field [29,33]. These observations can 
be explained by stimulus-invariant changes occurring in 
higher decision-making cortical areas that augment those in 
lower visual areas. 

In our model of learning- induced plasticity of VI, we did not 
consider if/how changes in VI could improve perception. In order 
to understand how perceptual learning could be generalized to 
other stimuli and locations, we would need to expand the model to 
include read-out decision-making circuits. However, modeling 
plasticity in the higher cortical areas, in addition to the changes in 
VI, is beyond the scope of this study, which is focused on changes 
that have been observed in VI. 



Conclusion 

Our study suggests that two previously reported inconsistent 
experimental results [7-8] may represent recordings from two 
different populations of VI neurons; One population of neurons 
that prefer the trained stimulus, and another that prefer different 
stimuli but still respond to the trained stimulus. 

Our model also predicts that learning results in strengthening 
feedback inputs into VI, which cause shunting inhibition in VI 
neurons. This novel mechanism of learning can be tested by 
measuring input resistances in relevant VI pyramidal cells during 
in vivo stimulation before and after learning. 

Despite of limited nature of learning-induced changes in V 1 , 
these changes could contribute to other learning-dependent 
changes in visual cortex, because of the hierarchical organization 
of the visual system [34—35]. Indeed, a learning-induced sharp- 
ening of orientation tuning curves was observed in an extrastriate 
visual area V4 [36] , which could be at least partially due to more 
sharply-tuned inputs from VI. 
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Figure 5. Stimulus specificity of tuning curve changes. A. Training the model with stimulus A strengthened the feedforward and feedback 
connections between stimulus-specific populations of VI and V2 neurons. B. A test stimulus B (green) activated population of neurons in VI and V2 
that included some neurons with the modified synapses. The red circle shows where the trained stimulus A was presented to the model. C. Another 
test stimulus C (purple) activated populations of VI and V2 neurons that did not include neurons with modified synapses. D. Before training, the 
strength of the feedback inputs was equal for any applied stimulus (grey dashed line). After training, the strength of feedback inputs depended on 
the applied stimulus (black solid line). E. Stimulus specificity of changes to tuning curves. Peaks of tuning curves did not change, but amplitudes and 
widths of the tuning curves changed depending on the distance from the trained orientation. 
doi:1 0.1 371 /journal.pcbi.1 003770.g005 



Feedback projections are ubiquitous in cortex [37-38], but 
despite recent progress [39-42], Kttle is known about functional 
roles of feedback connections in cortical plasticity. Our work 
presents a new hypothesis for the role of the feedback connections 
in cortical plasticity. We propose that top-down inputs are plastic 
and involved in perceptual learning. This prediction can be tested 
by blocking synaptic plasticity in extrastriate visual cortical areas 
and observing whether this would diminish learning-induced 
improvement of behavior. 

Models 

Spil<ing neuron networl< model of visual cortex 

We explore the feedback hypothesis in a spiking neural network 
model of the visual cortex (Fig. 2). The model consisted of two 
layers. Layer one represented the primary visual cortex VI and 
layer two represented a higher cortical area V2. 

In the biological visual cortex, feedforward inputs from V 1 drive 
activity in V2 neurons, but feedback inputs from V2 neurons back 
to VI do not drive but modulate activity in VI [43^4], thus 
avoiding a strong positive feedback loop that could lead to 
uncontrolled epileptic-like oscillations [45]. In the model, VI 
neurons activated V2 neurons through the convergent feedforward 
connections (Fig. 2). The diffuse feedback inputs from V2 did not 



drive VI neurons, but altered their response functions through 
shunting inhibition mechanism [1 1-12] . Shunting inhibition in the 
model was due to the balanced excitatory and inhibitory to the 
excitatory VI neurons, which modulated neural responses in VI 
neurons by subtracting neural responses to driving stimuli. 

VI layer. In the primary visual cortex, orientation selectivity of 
layer 2/3 neurons are formed by orientation- tuned feedforward 
inputs from layer 4 neurons, and by recurrent interactions among 
layer 2/3 neurons [15]. Orientation selectivity of layer 4 simple cells 
in turn is formed by feedforward inputs from arranged LGN cells 
[46]. In this model, we did not model LGN and layer 4 neurons. 
Layer 2/3 neurons receive a bell-shaped input, which represented 
orientation tuned inputs from layer 4 cells. The neurons in VI layer 
were organized with preferred orientations evenly spaced along one- 
dimensional network. Stimuli of different orientations were modeled 
by sliifting the center of the input along one-dimensional VI network. 

VI layer contained 100 excitatory and 25 inhibitory neurons, 
which represented layer 2/3 of VI. Both the excitatory and inhibitory 
neurons were simulated using Hodgkin-Huxley type neuron model 
with sets of parameters resulting in regular spiking activity [47]. 

Excitatory neurons. 



c 



dt 
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Figure 6. Effect of feedbacks on population response in V1. A. Population response in V1 to a stimulus before training, with (black) and 
without (grey) recurrent connections among VI neurons. B. Input-output function in VI neurons, with (black) and without (grey) recurrent 
connections among VI neurons. C. Population response in VI to the trained stimulus (preferred stimulus for neuron A) before and after training. D. 
Responses of neurons A and Sto the trained stimulus (preferred for neuron A) changed due to feedback. Black line indicates input-output relation in 
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VI. Blue arrows indicate the strength of feedback inputs. Black and grey filled circles represent neural response of neurons A and 5 to the trained 
stimulus before and after training. E. Population response in VI to a novel stimulus (preferred for neuron Bt, before and after training with the 
stimulus preferred for neuron A. F. Responses of neurons A and Bxo the novel stimulus (preferred for neuron B) changed due to feedback. Black line 
indicates input-output relation in VI. Blue arrows indicate feedback inputs. Black and grey filled circles represent neural responses to stimulus 
preferred for neuron B before and after training. 
doi:1 0.1 371/journal.pcbi.l 003770.g006 
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potassium current. 

Model parameters were chosen to fit the spike properties of 



= 0.05, gNa = 

Vn, = 55, 



100, 



gK- 



Vk 



40, 
80, 



regular spiking pyramidal neurons: gL, 
gc. = 0.9, gAHP = 0.05, Vl=-65 

Vca= 120, TCa= 100. 

Inhibitory neurons. The membrane potential of the inhib 
itory neurons was indicated as follows 



a„(F) = 



0.1-(F-l-30) 

1 _g-0.1(F+30) 



dh 

dt 



j6„(F) = 4-e-<''+«>/"* 



a,(F) = 0.07-e-(^+^V20 



hiV)- 



1 

1 _g-0.01(K+14) 



dn 
dt 



= 3-(a„(F)-(l-«)-iS„(F)-«) 



0.01-(F-H34) 

l_e-0.1(F+34) 



iS„(F) = 0.125-e-(^+^'/80 



lAHP=gAHP' 



[Ca 



l + [Ca] 



(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



d[Ca] _ -Omi-gCaiVe-Vca) 



dt 



l+e-(F«+25)/2.5 



■[Ca]/xca (13) 



where Vg is membrane potential of an excitatory neuron, m, h, 
n are the gating variables of voltage-gated sodium and 
potassium channels, gj^ is maximum conductance of the leak 
current, gi<fa is maximum conductance of sodium current, gx is 
maximum conductance of potassium current, Vl is reversal 
potential of the leak current, Fv„ is reversal potential for 
sodium currents, Vk is reversal potential for potassium currents, 
[CaJ is the intracellular concentration of free calcium ions, and 



dVj 

^ dt ~ ^^^^^ ^mem ^syni ( 1 4) 

where F, is membrane potential of an inhibitory neuron, and m, h, 
n are the gating variables for the ionic currents. The dynamics of 
the gating variables m, h, n were the same as for the excitatory 
neuron, the only differences between excitatory' and inhibitory 
neurons were, first, that the inhibitor)' neurons lacked the Iahp 
current and, second, in the inhibitory intc-rneurons membrane, 
time constants were shorter. Since the membrane time scale tmem 
is related to the membrane leak conductance gL as Tn,em~ l/gL> we 
have chosen in the inhibitor)- cells gii = 0.1 and in the excitatory 
cells gL,. = 0.05. Thus, the dynamics of the inhibitory cells was 
faster compared to the dynamics of the excitatory cells. 

Synaptic currents. Synaptic currents to an excitatory 
neuron were then calculated according to 



(15) 



Synaptic current to an inhibitory neuron: 



(16) 



where Vg is membrane potential of an excitatory neuron, F; is 
membrane potential of an inhibitory neuron, Sg is synaptic gating 
variables of AMPA currents, Si is synaptic gating variable of 
GABAa current, fg,fi are synaptic depression factors, and F,„ = 80, 

V,, = 0, V„ = 0. 

Synaptic models. Models for the AMPA receptor and 
GABA receptor were applied to simulate excitatory connections 
and inhibitory connections, respectively 



dse 
dt 



l-Se 



l+e-(F«+20)/4 



Se 1 / 'CaMPA 



I -Si 



dsi_ _ 

dt "ll+e-(n+20)/4 



— Si I /ZGABAa 



(17) 



(18) 



where Sg is a .synaptic gating variable of AMPA current, and s, is a 
synaptic gating variable of GABAa current, Xampa = 2 and 

GABAa — 2- 
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V2 layer. The V2 layer contained 100 excitatory neurons. To 
model the excitatory neurons in V2, we used Hodgkin-Huxley 
equations with the same parameters as those used for layer one 
excitatory' neurons. 

Connections within layers. The excitatory and inhibitory 
neurons in VI were connected through recurrent connections. 
Weights of the recurrent connections were distance-dependent and 
followed Gaussian distributions with s.d. = 10 neurons. The ratio 
between excitatory and inhibitory weights were chosen to balance 
recurrent excitation and inhibition in VI neurons. The recurrent 
connections did not change with learning. There were no recurrent 
connections among V2 neurons. 

Connections between layers. The excitatory neurons in VI 
layer projected to the excitatory neurons in V2 layer through 
convergent feedforward connections. Probability of feedforward 
connections depended on corresponding locations of neurons in VI 
and V2, and were randomly selected from Gaussian distributions 
with s.d. = 5 neurons, which result(;d in a compact stimulus 
representation in V2 layer. The feedback connections from excitatory 
V2 neurons projected uniformly to all (excitatory and inhibitory) 
neurons in VI. 

Learning. We simulated stimulus learning in the model by 

repeatedly presenting a stimulus with a particular orientation. 
Training the model strengthened excitatory feedforward connec- 
tions from neurons representing the trained stimulus in VI to 
neurons representing the trained stimulus in V2. Strong 
feedforward connections increased activity in V2 neurons, which 
in turn triggered plasticity of feedback connections. 

Plasticity of feedforward connections. Feedforward con- 
nections from the excitatory VI neurons to V2 neurons were 
adjusted according to a Hebbian learning rule. In the learning 
rule, we used intracellular calcium concentration as an indicator of 
a neural activity [48] 



where Xnebh =100 msec is time scale of Hebbian learning, W^^ is 
the strength of a feedforward connection from an excitatory VI 
neuron j to V2 neuron i, CaJ^ is intracellular calcium 
concentration in the excitatory VI neuron i, Ca^^ is intracellular 
calcium concentration in the excitatory V2 neuron j. 

Plasticity of feedback connections. The feedback projec- 
tions allowed V2 neurons to decrease activity in VI when neurons in 
V2 were activated strongly by feedforward inputs from VI . In order 
to achieve this modulatory effect in the model, the strength of the 
feedback connections were adjusted according to a plasticity rule 



-Can 



(20) 



•J K2 ^ VI 

T^Hebb ■ — — = Ca,. -COj 



where Tfb = 500 msec is a time scale of feedback plasticity, Gf* is 
the strength of the feedback connection from V2 neuron i to the 
excitatory and inhibitory VI neurons, Caf^ is an intracellular 
calcium concentration in the excitatory V2 neuron i, IS a 

homeostatic value of the intracellular calcium concentration in the 
excitatory V2 nc'urons, which was calculated as a time averaged 
value of the intracellular calcium concentration in the excitatory V2 
neurons. 
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